Counting Lattice Points by means of the Residue Theorem 1 
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Abstract. We use the residue theorem to derive an expression for the number of lattice 
points in a dilated n-dimensional tetrahedron with vertices at lattice points on each co- 
ordinate axis and the origin. This expression is known as the Ehrhart polynomial. We 
show that it is a polynomial in t, where t is the integral dilation parameter. We prove the 
Ehrhart-Macdonald reciprocity law for these tetrahedra, relating the Ehrhart polynomials 
of the interior and the closure of the tetrahedra. To illustrate our method, we compute 
the Ehrhart coefficient for codimension 2. Finally, we show how our ideas can be used to 
compute the Ehrhart polynomial for an arbitrary convex lattice polytope. 



1. Introduction 

Let Z n C R n be the n-dimensional integer lattice, and 

V=\ (xi, . . . , x n ) € R n : V — < 1 and all x k > 1 
I ti ak J 

the n-dimensional open tetrahedron with vertices (0, . . . ,0), (a%,Q, . . . ,0), 
(0, d2, 0, . . . , 0) , . . . , (0, . . . , 0, a n ). where a%, . . . , a n are positive integers. 
For t € N, denote by L(V, t) and L(V, t) the number of lattice points in the 
dilated polytope tV and its closure, respectively. 

Ehrhart ( |Ehj ) proved (actually for a general lattice polytope) that both 
L(V, t) and L(V, t) are polynomials in t of degree n. Moreover, he deter- 
mined the two leading coefficients and the constant. The leading coefficient 
is Vol('P), and the second coefficient is 1/2 Vol(dV), which is half the sur- 
face area of V normalized with respect to the sublattice on each face of V . 
The constant coefficient equals x('P), the Euler characteristic of V. The 
other coefficients of L(V,t) and L(V,t) are not as easily accessible. In fact, 
a method of computing these coefficients was unknown until quite recently 
( |Baj . |Br-Vej . |Di-Ro| . |Ka-Khj . |Kh-Puj ) . In this paper, we present an 
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elementary method for computing the Ehrhart polynomials of V and V us- 
ing the residue theorem. We verify the Ehrhart-Macdonald reciprocity law 
for these ra-dimensional tetrahedra. To illustrate our method, we compute 
the first nontrivial coefficient, c ra _2, of the Ehrhart polynomial. Finally, we 
show how our ideas can be used to compute the Ehrhart polynomial for an 
arbitrary convex lattice polytope. 



2. The main idea 

Let's start with L(V, t); that is, we consider the closure of our dilated tetra- 
hedron tV. We introduce the notation 

A := a\ ■ ■ ■ a n , Afr := ai ■ ■ ■ a,k ■ ■ ■ a n , 

where a& means we omit the factor a&. Then we can write 



L(P, t) = card I (mi, . . . , m n ) € Z n : — < t and all m k > I 
I fc=i ak J 

= ca r d{ (mi ,...,,„„, m)eZ »«: %£?* + m = tA }■ 

We can interpret L(V,t) as the Taylor coefficient of z tA for the function 
(1 + Z M + Z 2A, + ^^ l + z M + Z 2A 2 +...)... 
• (l + z An + z 2An + . . .) (l + z + z 2 + . . .) 



11 11 



1 - z A ^ 1 - z A ? \-z A ^\-z 



Equivalently, 



_ / z~ tA ~ l \ 

t] = Res U-**)(i-**)---(i-^)(i-*) '* = v • 



To reduce the number of poles, it is convenient to change this function 
slightly; this residue is clearly equal to 

/ z -tA _i \ 

ReS V(l - z A ^) (1 - z A *) •••(!- z^) (l-z)z' z = ) +1 - 
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If this expression counts the number of lattice points in tV, then all we have 
to do is compute the other residues of 

z -tA _ i 

:= (1 - z A ^) (1 - z A *) ■■■{!- z A -) (l-z)z 
and use the residue theorem for the sphere C U {oo}. In this notation, 

L(V,t) = Res (/_*(*), z = 0) + l ■ (1) 
The only poles of /_< are at 0, 1 and the roots of unity in 

tt := \z e C\ {1} : z^i = 1,1 < k < j < n\ . 



Note that Res(/_t,z = oo) = 0, so that the residue theorem gives us the 
first half of our main result: 

Theorem 1(a). 

L(P, t) = l- Res (f- t (z), z = 1) - R es (/-*(*). « = V 

xen 

Remarks. 1. The residue at z=l can be calculated easily: 

Res (f-t(z), z = l) = Res (e z f_ t (e z ), z = 0) 
/ e -tAz _ 1 

= ReS V (1 - e A ^ z ) (1 - e A * z ) • • • (1 - e A « 2 ) (1 - e z ) ' ^ = ° 

To facilitate the computation in higher dimensions, one can use mathe- 
matics software such as Maple or Mathematica. It is easy to see that 
Kes(f-t(z), z = 1) is a polynomial in t whose coefficients are rational ex- 
pressions in ai, . . . , a n . 

2. The residues at the roots of unity in 0, are in general not as easy to 
compute. They give rise to Dedekind-like sums and their higher dimensional 
analogues, as we will illustrate in section 4. There is, however, one feature 
we can read off from these residues immediately, the dependency on the 
dilation parameter t: 

Corollary 2 (Ehrhart). L(V,t) is a polynomial in t. 

With Corollary 3 below, this will also imply that L(V,t) is a polynomial. 
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Proof. Let A € f2 be a l?'th root of unity, where B is the product of 
some of the a k - Now express z~ tA in terms of its power series about z = 
A. The coefficients of this power series involve various derivatives of z~ tA , 
evaluated at z = A. Here we can introduce a change of variable: z = vjb = 
exp (-^ log w) , where we choose a suitable branch of the logarithm such that 
exp (-jjjlog(l)) = A. The terms depending on t in the power series of z~ iA 
consist therefore of derivatives of the function z~ tA l B , evaluated at z = 1. 
From this it is easy to see that the coefficients of the power series of z~ tA 
are polynomials in t. The fact that L(V, t) is simply the sum of all these 
residues, finally, gives the statement. □ 

For the computation of L("P, t) (the number of lattice points in the interior 
of our tetrahedron tV), we similarly write 

L(V, t) = card I (m u ... , m n ) G Z n : V — < t and all m k > 1 
I ti ak J 

= card((m 1 ,...,m„,m)G^+ 1 : ^U^A k + m = tA \ 
Now we can interpret L(V, t) as the Taylor coefficient of z for the function 

(^i + + ) ( Z A 2 + Z 2A 2 + _)... + Z 2A„ + _ ) , z + z 2 + ) 

z Al Z M z An z 



1 - z A ^ 1 - z A * 1 - Z A " 1 - z ' 
or equivalently as 

_ / Z Al Z^ 2 2 

Res 2- — ■ ■ ■ — - z \ z = 

\ 1 — z Al 1 — z 712 1 — z An 1 — z 

( z Al z A2 z An z z~ tA - 1 

= Res T^z^ " ' 1 - z A ™ T^~z z ,Z = ° 

To be able to use the residue theorem, this time we have to consider the 
function 

11 11 z tA — 1 



z Al — 1 z A ' 2 — 1 z An — 1 z — I z 
so that 

L(P,t) = (-l) n Res(f t (z),z = oc) . (2) 
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The finite poles of ft are at (with residue -1), 1, and the roots of unity in 
f2 as before. This gives us 

Theorem 1(b). 

L(V, t) = (-If ( 1 - Res (f t (z), z = 1) - Res (f t (z),z = A) 

As an immediate consequence we get the remarkable 

Corollary 3 (Ehrhart-Macdonald Reciprocity Law). 

L(V,-t) = (-l) n L(V,t) . 

This result was conjectured (again, for a general lattice polytope) by Ehrhart 
(jEEl), and later proved by Macdonald (|Maj), McMullen ([Mel), and Stanley 



3. The Ehrhart coefficients 

With a small modification of ft(z), we can actually derive a formula for each 
coefficient of the Ehrhart polynomial 

L(V,t) = c n t n + ... + c . 

Consider the function 
9k{z) ■-- 



tA -\* 



(1 - z A i) (1 - z A *) ■■■{!- z A ") (l-z)z 

E* =0 Q)z-^)(-iy 



(1 - z A ^) (1 - z A *) ■ • • (1 - z A ™) (1 - z) z 
If we insert — X]j=o Cj) ( — ^¥ = m the numerator, this becomes 
1 'k\, 4 , z -t(k-j)A_ 1 



9k{z) = 

j=0 



^ G) ( (1 - Z A ^) (1 - Z A *) ■ ■ ■ (1 - Z A ") (l-z)z 

E0(-i) j /- t H)(^ 



5 



Recall that (JTJ gave us L(V, t) = Res (f-t{z), z = 0) + 1. Using this relation, 
we obtain 



Res (efc(z), * = 0) = J2 ( • ) ^s (/-t^*), * = 0) 

3=0 

k — 1 

=e(-) ( - i)j ( l ^ (a; -^)- 1 ) 

i=o v/ 

E(*)(-l^^(*-^) + ("I)"- 



3=0 

We claim that this polynomial has no terms with exponent smaller than k: 
Lemma 4. Suppose L(V, t) = c n t n + . . . + cq. Then for 1 < k < n 

n 

Res (g k (z), z = 0) = k\ ^ S(m, k) c m t m , (3) 

m=k 

where S(m, k) denotes the Stirling number of the second kind. 



Proof. Suppose 

fc-i 



//A n 

E • ( fe - = E 6 *.™* m > 



(4) 



j=0 yj/ m=0 

so that for m > 

6 *.™ = E (*) c- 1 )^* - i) m = c - E Q) (- r 



fc-i -m 



The Stirling number of the second kind S^m, A;) is the number of partitions 
of an m-set into k blocks. The reason we are interested in these numbers is 
the identity ([Si]) 



3=0 KJ/ 



-i) k - j r, 



so that bk, m = Cm k\ S(m, k) for m > 0. Some of the elementary properties 
of S(m,k) are (|S2) 

S(m, k) = if k > m (5) 
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5(m,l) = l (6) 
S(m, m) = 1 (7) 
S(m, k) =k S{n -l,k) + S(n - 1, fc - 1) . 



By (JSJ), we conclude that bk t m = for 1 < m < k. The constant term in 
is 

i— n \J/ 



j=0 



Since Co = 1 for our tetrahedron (in fact, cq = 1 for any convex polytope 
[EE]), © follows. □ 



The other poles of gk are at 1 and the roots of unity in 

A 



n k := iz G C \ {1} : z a ^- a ^+i = 1, 1 < j x < j 2 < . . . < j fe+1 < n | . 

Note that as /c gets larger, fi^ gets smaller. That is, we have fewer residues 
to consider. This is consistent with the notion that the computational com- 
plexity increases with each additional coefficient, that is, the computation 
of Cfc is more complicated than that of Ck+i- Using the residue theorem, we 
can rewrite © as 

Theorem 5. Suppose L(V, t) = c n t n + ...-)- cq. Then for 1 < k < n 

S(m, k) c m t m = ^ j Res (g k (z), z = 1) + Res (g k (z), z = A) J . 

m=k ' \ Aen fc J 



Remarks. 1. For k = 1, we get with Q a reformulation of Theorem 1(a). 

2. The coefficients of L(V,t) are the same as those of L(V,t), up to the 
sign: By Corollary 3, L(V,t) = c n t n - c n ^t n - x + ... + (-l) n c . 

3. Res(gk(z), z = 1) can be computed as easily as before, the slightly more 
difficult task is to get the residues at the roots of unity (see also remark 2 
following Theorem 1(a)). However, with increasing k, we have to consider 
fewer of them, so that there is less to calculate. If we want to compute the 
Ehrhart coefficient c m , we only have to consider the roots of unity in O m . 
We can make this more precise: With (j7|, we obtain 
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Corollary 6. For m > 0, c m is the coefficient of t m of 



— I Res (g m (z), z = 1) + ^ Res (g m (z),z = A) 



4. An Example 

As an application, we will compute the first nontrivial Ehrhart coefficient 
c n _2 for the ra-dimensional tetrahedron V (n > 3) under the additional 
assumption that ai,...,a n are pairwise relatively prime integers > 2. This 
case was first explored by Pommersheim ( \Po\ )- 

Theorem 7. Under the above assumptions, 

Cn-2 = ^ - 2 y (c n - s{Ai,a 1 ) - ... - s(A n , a n )) , 
where s(a,b) denotes the Dedekind sum, and 

1 , 1 / 1 M K 

C n := 7 (n + Ai, 2 + . . . + A n _i, n + — T + — + ... + — 

Here A,-^ denotes a\ ■ ■ ■ a,j ■ ■ ■ at ■ ■ ■ a n - 
Proof. We have to consider 



9n-2{z) 



(1 - (1 - Z^) •••(!- *^») (I- z)z 



Because oi, . . . ,a n are pairwise relatively prime, g n -2 has simple poles at all 
the ai, . . . , a n th roots of unity. Let A ai = 1 ^ A. Then 

Res (g n -2(z),z = A) = 

1 / (z- tA -l) n - 2 \ 
Res - ^ J——,z = \). 



(1 — X Al ) (1 — A) A Ul - z A 2) ... (l - ' 

Similar to the methods used in the first chapter to arrive at Corollary 2, we 
make a change of variables z = w l / ai = exp ( — log w ) , where we choose a 
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suitable branch of the logarithm such that exp ^log(l)^ = A. We thus 
obtain 

Res (g n - 2 {z),z = A) = 

Res T, 77 R"T) w = 1 > 



" (l-A A i)(l-A)Aai Wl-w; B 2)---(l-w; B «) ; 
where B := 02 ■ ■ ■ a n , := 02 ■ • • • • ■ a n . We claim that 
R es f ( z _ \ _ 2 

KeS ^(l-^)---(l-^)' Z " J " * • 

To prove this, first note that 

[z- tB - l) n ~ 2 = (-tB) n -\z - If" 2 + 0((z- If' 1 ) . 
Now for m G N, 

« / 1 \ 1 «-! 1 

Res , z = 1 = lim 



l-z m y z-*il-z m m 

Putting all of this together, we obtain 

^-rv.^ ( -^ 2 



(l-^)---(l-^)' y (_ jB2 )...(_ jBn ) 

= -f- 



+n—2„n—2 „n—2 
i Uo • • • u„ ,n—2 



Q>2 ' ' ' ® n 

as desired. Therefore 

_ f n-2 



Res (g n - 2 (z),z = A) 



ai (l-A^)(l-A) 
Adding up all the ai'th roots of unity 7^ 1, we get 



_-f-n—2 j 

£ Ba>Gfc- 2 (*),* = A) = — £ (1 _ A A 1)(1 _ A) 



A a i=l^A A a i=l^A 
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where £ is a primitive cti'th root of unity. This finite sum is practically a 
Dedekind sum: 



«i jfe (i-e fcj4i )(i-e fc ) 4«i ^ v i-e M ) v i-ev 

1 . , . i <l ^-^ 1 / 7T &t4i 7rA; \ 

= - — (ai - 1) - - — > cot h cot — 

4ai 4ai ^ \ a i a i / 

1 a ^-^ 7rA;yli 7rfc 

— - — > cot cot — 

4ai f-f ai ai 

fc=i 

= I-4^- 5( ^' ai) - 

The imaginary terms disappear here, since the sum on the left hand side 
and s(Ai,ai) are rational: Both are elements of the cyclotomic field of ai'th 
roots of unity, and invariant under all Galois transformations of this field. 

Hence we obtain 

Res (g n - 2 (z),z = A) = -f 1 ' 2 ( ± - -L - 8 (A U ai ) ) . 

A«i=l^A V4 4fll 7 

We get similar expressions for the residues at the other roots of unity, so 
that Corollary 6 gives us for n > 3 

1 (n 1/1 1 

-s(i4i,ai) - ... -s(A n ,a n )\ (8) 

where C is the coefficient of t n ~ 2 of Res (g n - 2 {z), z = 1). We can actually 
obtain a closed form for C: As before, 

Res (g n - 2 (z), z = I) = Res {e z g n - 2 (e z ), z = 0) 

/ (e^'-l)"" 2 \ 

" ReS ^ (1 - e A ^) (1 - e A **) • • • (1 - e A ^) (1 - e z ) ' * ~ °) ' 

Now with 

( e -t^ _ i) n " 2 = (_t^) n - 2 + O ((te) n_1 ) 
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the coefficient of t n 2 of Res (g n -2(z), z = 1) turns out to be 





4 \A 2 ■ ■ ■ An Ax - • • A n _i 



(~l) n+ M n 

Ai • • • A n _i 



+ 



+ 




(-1) 



n-l 



) 



Substituting this into (JSJ) yields the statement. 



□ 



The other Ehrhart-coefficients for this tetrahedron can be derived in a sim- 
ilar fashion, although the computation gets more and more complicated, as 
noted in section 3. 



5. General lattice polytopes 

Any convex lattice polytope (that is, a convex polytope whose vertices are 
on the lattice Z n ) can be described by a finite number of inequalities over 
the integers. In other words, a convex lattice polytope V is an intersection 
of finitely many half-spaces. Translation does not change the lattice point 
count, so we can assume that the points in the polytope have positive co- 
ordinates and apply the ideas of the previous sections to V . Suppose the 
closure of the dilated polytope tV is given by the n + q inequalities 



X\ , . . . , x n ^ 




0) 
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with cijk E Z. Then by introducing a similar notation as before, 

Pk ■ — O'kX ' ' ' O'km -"-jk ■ — CLjl ' ' ' Qjk ' ' ' ^jm 

we can define a matrix 



M = ( A jk 



3=1-9 



Let Cj the j'th column, and R k the fc'th row of M. Then we can rewrite 
the last q inequalities of © as 

Ri-x< tPi 

(10) 

where x = (x\, . . . , x n ) and • denotes the usual scalar product. Now consider 
the function 

v - t p x -\ _ v -tP q -\ 

f( z ) = f(zi,---,z q ) 



{ l- z Ci ) ... { i_ z c q){ i_ Zx) ...( 1 _ Zq y 

Here we used the standard multinomial notation z v := z^ 1 ■ ■ ■ z"q . We will 
integrate / with respect to each variable over a circle with small radius: 

f(z 1} ...,Zq) dZg-'-dZ! (11) 

\Zg\=€ q 

Here, < eo, . . . , e q < 1 are chosen such that we can expand all the 1 _ 1 z c k 
into power series about 0. To ensure the existence of eo, • • • , e q , we may, if 
necessary, add an additional inequality x\ + . . . + x n < tPo for a suitable 
large Pq. This is always possible, since V is bounded. 

Since the integral over one variable will give the respective residue at 0, we 
can integrate with respect to one variable at a time. When / is expanded 
into its Laurent series about 0, each term has the form 

m-Ri+n-tPi-l m-R q +r q -tP q -l 

z i ' ' ' z g ) 

where m := (mj, . . . , m n ), and mi, . . . , m n , r±, . . . , r q are nonnegative inte- 
gers. Thus, in the integral Qll|) . this term will give a contribution precisely 
if m satisfies the inequalities (|10[). In other words, we have proved 
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Theorem 8. 

L(P, t) = -77T-7T- / ■•• / f(z 1 ,...,z q )dz q ---dz 1 . 
(2vrz)9 J\ Zl \ =tl J\z q \=e q 

On the other hand, we can compute this integral by computing the residues 
/ as before, for each variable at a time. The (finite) sum over these 
residues will be equal, by Theorem 10, to the number of lattice points in the 
closure of our polytope tV. We will explore this method and its results in a 
future paper. 



6. Further remarks 

1. The ideas in this paper can also be used in other ways. First, note that 
by omitting the z ±tA — 1 factor, we essentially get trivial residues at oo and 
0. In view of Corollary 6, this corresponds to obtaining an expression for 
the constant Ehrhart coefficient and setting it equal to 1 ( |Bej ). This means 
we are not getting a lattice point count; however, the other factors will still 
give us Dedekind-like sums, as they appear in the formulae. By using the 
residue theorem as before, we thus now get relations between these terms, 
similar in spirit to the various reciprocity laws for classical ( |R,a-Gr| ) and 
higher dimensional Dedekind sums ( |Zaj ) . 

2. Finally, we can extend these methods to n-dimensional rational polytopes 
(whose vertices have now rational coordinates). For example, the number 
of lattice points in the dilated rational tetrahedron 

1 

(xi ,x n ) 6 R n : x kik < t an d all > > 
k=l J 

can be computed with similar methods. The details will appear in a future 
paper ( |Be-Di-Ro] ) . 
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